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ABSTRACT 

We monitored Eta Carinae with HST WFPC2 and Gemini GMOS throughout the 
2009 spectroscopic ev ent, which was expected to differ from its predecessor in 2003 



([Davidson et al.l 12005). Here we report major observed differences between events, and 
their implications. Some of these results were quite unexpected. (1) The UV brightness 
minimum was much deeper in 2009. This suggests that physical conditions in the early 
stages of an event depend on different parameters than the "normal" inter-event wind. 
Extra mass ejection from the primary star is one possible cause. (2) The expected 
He II A4687 brightness maximum was followed several weeks later by another. We 
explain why this fact, and the timing of the A4687 maxima, strongly support a "shock 
breakup" hypothesis for X-ray and A4687 behavior as proposed 5-10 years ago. (3) We 
observed a polar view of the star via light reflected by dust in the Homunculus nebula. 
Surprisingly, at that location the variations of emission-line brightness and Doppler 
velocities closely resembled a direct view of the star; which should not have been true 
for any phenomena related to the orbit. This result casts very serious doubt on all the 
proposed velocity interpretations that depend on the secondary star's orbital motion. 
(4) Latitude-dependent variations of H I, He I and Fe II features reveal aspects of wind 
behavior during the event. In addition, we discuss implications of the observations for 
several crucial unsolved problems. 
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Introduction 



Eta Carinae is exhibit A for episodic mass loss near the top of the H-R diagram, for the physics 
of giant eruptions ("supernova impostors") and subsequent recovery, the behavior of outflows above 
the Eddington limit, polar winds, several exotic nebular processes, etc. Most of these topics remain 
poorly understood, but there is no reason to think that rj Car structurally differs from other 
extremely massive stars. A likely companion object affects the phenomena, but does not by itself 
constitute an "explanation." 

Beginning in the mid-1940s rj Car began to exhibit occasional spectroscopic changes that we 
now recognize as a 5.5-year spectroscopic/photometric cycle. Occasio nally its high-excitation He I, 



[Ne III], [Fe III] emission lines disappear for a few weeks or months (Gaviola 



19841) while other changes also occ ur, specifically in the X - ray (e.g., 



1999 



and infrared flux (e.g., IWhitelock et al 



199 



4; 



1953: 



Corcoran et al 



Fe ast et al 

events" recur with a period close to 2023 davs (Ipaminerilll9*9t] : 



2001 



Zanella et al 



1997; 



Ishibashi et al 



These "spectroscopic 



1994 



Damineli et al 



— .Wh itclock et a 

1999 . l200(ilwiutelock et akhooi IMartin et al.ll2006al : lDamineh et all2008bl 

2010 ). They have been attributed to (1) eclipses of a hot secondary star by the primary win d 



Fernandez-Lajus et al 



(Damineli et al. 



1997 



Ishibashi et al. 1999b; Stevens fc Pittard 



1999 



Pittard Corcoran 



(2) d i sturbances in the primary wind trigger ed by a companion star near periastron ( Davidsonl 



1999; Smith et al 



2003 



19841; Davidson et al 



Martin et al 



2000; 



Smith et al 



200 



2006a); ( 3 ) a thermal/rota tional recovery cycle (| Zanella et al 



2) 



1997 



20031 ; iDavidsonl 120051); or (4) a breakup /collapse of the 



wind-wind collision structur e due to known shock instabilities ( Davidsonll2002 ; Soker|[2003 ; Martin et al 
2006al ; ISoker fc Beharl 120061 ) . These ideas are not mutually exclusive. Either (2) or (3) would be 
significant for massive-star physics in general, because they may require an undiagnosed instability 
near the Eddington limit. 

Observations of the 2003.5 event appeared to favor possibility 2 and especially 4, but did not 
rule out number 3, and the likely geometry suggests that an eclipse probably occurred with lesser 
conseque nces. Photometric behavior, a chaotic X-ray b ehavior, and an unpredicted He II A4687 
outburst (jSteiner &; Daminefil 120041 ; iMartin et al.ll2006al ) were especially significant. 



Mean while the longer-term b e havior changed dramatically . The central star br i ghten ed rapidly 
after 1998 ([Davidson et al.lll999al ; IMartin & Koppelmanll2004l ; IMartin et al]|2006bl. boioh ■ an d ma- 
jor spectral features differed between the 1998.0 and 2003.5 events (jDavidson et al,ll2005l ). De- 
struction and/or lessened formation of dust played a role, but the root cause must involve a secular 
change in the UV flux or the wind density or both. Thus, from the viewpoint of 2007-2008, ob- 
servations of the expected 2009.0 event merited a high priority for comparisons with 2003.5 and 
1998.0. Unfortunately the Space Telescope Imaging Spectrograph on the Hubble Space Telescope 
(HST STIS) had failed in 2004, and thus was not available to separate the star from ejecta (see 
below). Later, after the event was observed as discussed in this paper, STIS became operational 
again and proved that emission lines from r/ Car's wind had greatly weakened since 2004; the mass- 
loss rate had probably decreased by a factor of the order of 2 or 3 in a time span of several years 
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(|Mehner et al.ll2010bl : see also iKashi & Sokerll2009al : ICorcoran et all 120101 ) . Therefore the 2009.0 
spectroscopic event occurred in physical circumstances appreciably different from its predecessors. 
This should be helpful for deducing the nature of the event, on the same principle as varying 
parameters in an experiment. 

Ground-based spectra of r\ Car represent an unresolved mixture of the central star plus bright 
ejecta located at r ~ 072 to 2". Fortunately, though, the central star has brightened more than 
the nearby ejecta in the past decade, so ground-based observations of it have become relatively 

contacted than at earlier toesEl Thus, in 2007 we began to observe the star and several 
offset positions in the Homunculus nebula with GMOS on the Gemini-South telescope. Apart from 
the question of secular changes, we planned to measure some aspects of the spectroscopic event 
better than had been done in 2003 and earlier. A primary goal was to monitor the behavior of the 
peculiar He II A4687 line with more frequent observations during the 2009 event. We also observed 
the spectrum reflected by dust in the southeast (SE) Homunculus lobe, a "polar" view of the star. 

In this paper we discuss the resulting spectra and the lig ht curve of the central s tar. Some 
aspects of the 2009 event have been report ed by other authors: Richardson et al.1 (|2010i ) described 



the behavior of H a, IKashi &: Sokerl (|2009al ) and ICorcoran et al 



Groh et al 



(201 0j) c o mmen ted on the X-rays, 



(2010 1) reported hig h- velocity material. iFernandez-Laj us et al.1 (|2010l ) discussed photom- 



etry, and iTeodoro et al.1 (120111 ) have recently discussed the He II A4687 emission. Here our scope 
is broader. As we explain below, some of the differences compared to the 2003 event provide unex- 
pected new evidence for specific phenomena; UV photometry of the central star indicates changed 
conditions; reflected spectra showing the polar view appear inconsistent with published models 
of the velocities; and we strongly disagree with Teodoro et al. concerning the past observational 
record and theoretical interpretations. Our main new conclusions are that (1) differences between 
the 2003 and 2009 events give valuable and specific clues to the phenomena; (2) a "shock breakup" 
scenario proposed a decade ago now seems almost inescapable; (3) the observed Doppler velocities 
are far less straightforward than most authors have assumed; and (4) realistic quantitative models 
- as opposed to qualitative scenarios or idealized simulations - are badly needed but will be very 
difficult to achieve. We also comment on many other observational and theoretical factors. Some 
important ideas that were discussed in connection with the 2003 event remain valid. 

HST photometry and Gemini GMOS observations are described in the next section. The 
2008-2010 light curve of the central star and the unusually deep minimum during the 2009 event 
are discussed in Section 3. In Section 4 we discuss the peculiar He II emission during the event, 
its connection with the X-ray flux curve, and the "FOS4" polar view of He II emission compared 
with our direct line-of-sight view. The changing wind geometry during the 2009 event is discussed 
in Section 5. In the discussion section we summarize the results, emphasize the new information 
given by the 2009 event, and raise some outstanding questions. 



1 Here "central star" really means the opaque primary wind. The secondary star is also included, but its spectrum 
is too faint to be observable at accessible wavelengths. 
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Observations and Data Reduction 



2.1. HST Photometry with WFPC2 and STIS 



We have monitored the brightness changes of the central star in several ba nd-passes with pho- 
tometry from HST ACS/HRC a nd WFPC2 images and STIS spectra since 1998 (jlVIartin & Koppelman 
20041 ; iMartin et al.ll2006bl . |2010| ). During the 2009 event we monitored the brightness of the central 
star with the HST WFPC2 camera using the F255W and F336W filters. F255 W samples the 
NUV "Fe II forest" (jCassatella e t al.lll979l: l Altamore etal1ll986l: Viotti et al.lll989h. whic h greatly 



increases in opacity during a spectroscopic event (IDavidson et al 



1999b 



Gull et alJl200(l ). F336W 



includes the Balmer continuum augmented by various emission lines. These filters have been cali- 



brate d for direct comparisons with the HST ACS/HRC F250W and F330W filters (jSirianni et al 



2005). Our own careful checks have led us to conclude that comparisons in these filters across 
instruments are valid for r/ Car. The images were reduced using the standard STScI data reduction 
pipeline. Calibrated fluxes were measured in a / .'3 diameter weighted aperture following proce- 
dures described in our previous papers which, combined with the spatial resolution of the HST, 
minimize the influence from nearby bright ejecta. ACS-equivalent photometry was also synthesized 
from HST STIS data before mid-2004 and after mid-2009; spectra were extracted with a weighted 
parabolic cross dispersion profile similar to the virt ual aperture used to measure ACS/HRC images, 
convolved with the filter functions, and integrated (jMartin et al.ll2006bl ). The WFPC2 photometry 
from 2008.7 to 2009.3 and STIS synthetic photometry from 2009.6 to 2010.6 is listed in Table E 
Earlier data from ACS/HRC and WFPC2 images and STIS spectra can be found in the papers 
cited above. 



2.2. Gemini GMOS Observations 

To cover the 2009.0 event, we obtained ground-based slit spectroscopy of rj Car with the 
Gemini-South Multi-Object Spectrograph (GMOS) beginning in 2007 June through 2010 January. 
In most cases, we used the B1200 line grating at three tilt angles to cover the spectrum from A3700 
to A7500 A. The tf.'5-wide slit, oriented with a position angle of 160°, was placed at four different 
positions: on the star, two offset positions ±0 / .'75 relative to the star, and at a position known as 
"FOS4," 4'.' 5 from the star in the SE lobe of the Homunculus. (Operationally this was done b 
offsetting the slit —2" parallel to itself.) The star's polar spectrum is reflected by dust at FOS4J3 

Since most of this paper is concerned with the He II A4687 emission and features at nearby 



2 The name "FOS4" originated when it was a target for th e HST Faint Object Spectrograph in the 1990's 
( Davidson et al. 19951 : Humphreys et al. 19991 : Zethson et al. 1999h . Apart from small pointing differences, FOS4 is 
the same as the "Center of SE Lobe" in Figure 2c o f Hu mphreys et a l . and Figure 3c of Zethson et al. It was further 
discussed bv lSmith et~aH (|2003h . IWeis et al l l|2005h . and Istahl et~ZI l|2005h . The 1991-1997 HST FOS data may be 
worth re-examining today for comparison with STIS observations made in 1998-2010. 
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wavelengths, we concentrate on the blue spectra of the star and FOS4. Table [2] is a journal of those 
observations. True slit positions vary slightly from the four locations described above, due to a 
combination of blind offset errors and differential atmospheric refraction, which were comparable. 
No corrector was available, and our observing goals usually did not allow observations with a vertical 
slit orientation ( "the parallactic angle" ) . Therefore in each case we calculated the apparent position 
of the star as a function of wavelength, and used the offset slit position that was closest to the star 
for blue wavelengths. This procedure was adequate because, in effect, the slit positions overlapped 
due to the combination of slit width and seeing; we did not attempt to measure accurate absolute 
fluxes. For more information see Technical Memo 14 on the t] Car Treasury Program websitejfl 

We prepared two-dimensional spectrograms using the standard GMOS data reduction pipeline 
in the Gemini IRAF pa ckage, and extracted one-dimensional spectra via a routine developed earlier 



for use with HST STIS ( Martin et al.ll2006al ). At each wavelength our software integrates the counts 
along a line perpendicular to the dispersion, weighted by a mesa-shaped function centered on the 
local spectral trace. We used a mesa function with base-width = 11 pixels and top-width = 7 pixels, 
about Of! 8 and 0ff5 respectively. The seeing was roughly Of! 5-1'! 5, so each GMOS sp ectrum discussed 



represents a region about 1" across. The spectra were rectified using a LOESS fit (jCleveland 



1979 



Cleveland fe Devlinlll988l ). The pipeline wavelength calibration was improved using the interstellar 
Ca II K absorption line at A3935 A and the interstellar absorption line at A5782 A. The absolute 
wavelength scale was obtained with HST STIS spectra that have better wavelength calibrations. 
Throughout this paper we quote vacuum wavelengths and heliocentric Doppler velocities. 



2.3. Concerning Times and Phases 



When referring to "phase" in the 5.5-year cycle, we consistently use the Treasury Project 
system with period 2023.0 days, which has been employed since 2003 without any need for revision; 



see |http:/ /etacar.umn.edu/ [ comments at the end of Section 2 in iMehner et al.l l2010al . and an 
Appendix to this paper. In this paper, time within a spectroscopic event is denoted by t, such that 
t = at MJD 54860.0 (2009 Jan 29), MJD 52837.0 (2003 Jul 17), etc. Periastron mo st likely occurs 
within a few days of t = 0, but that was not the reason for the choice of zero point. iTeodoro et al. 
(|201ll ) arbitrarily use a different zero point that corresponds to t ~ — 18 d, but the rationale for 
the long-standing Treasury Project system is noted in the Appendix. 



3. Light Curve of the Star 2008 2011 

The light curve of the central star from HST data has two crucial advantages that have been 
overlooked in some discussions of r/ Car's photometric record. First, all ground-based photometry 



3 http:/ /etacar. umn.edu/treasury/techmemos/pdf/tmemo014.pdf 
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includes ejecta at radii r > (/.'15, practically unrelated to short-term variations of the central 
star. This contamination is remarkably strong, varies with time differently than the star, and 
is difficult to quantify. In recent years the ejecta accounted for 40- 70% of the br i ghtnes s in the 
best ground-based photometry at visual wavelengths; see Figure 3 in iMartin et al.l (|2006bl ). Since 
this relative fraction has lately been decreasing on a timescale of ~ 10 years (see the ref. just 
cited), future ground-based photometry may eventually become representative of the star but pre- 
2005 measurements were dominated by ejecta. The co ntamination was prob ably still important 
during the 2009 event, but the amount is unclear; see iMehner et al.l (|2010bl ) and the discussion 
below. Fortunately the high spatial resolution of HST allows us to sample the central object itself, 
r < / .'15, with little extraneous contamination - alb eit wi th less precision than high-quality ground- 



based photometry (cf. iFernandez-Laius et al. 



2009 



20101 and http:/ /etacar.fcaglp.unlp.edu.ar/). A 



second obvious and important advantage of HST is access to the UV region A < 300 nm, where 
the largest photometric changes occur during a spectroscopic event (see refs. in Section [2TT|> . 

Indeed the HST UV data illustrate an important point of this paper that is not obvious in 
the ground-based photometry: r\ Car's 2009.0 event differed from its 2003.5 predecessor. Figure 
DQ shows the light curve in the F250W and F330W filters from 1998 through 2010 based on HST 
ACS/HRC and WFPC2 images and STIS spectra. The two smaller panels expand the light curve 
around the time of the event and show the minima for the 2009.0 and 2003.5 events superimposed. 
The depth of the minimum in the 2009 event was about 1.1 mag in the F250W filter, which is 
sensitive to the Fe II forest, and 0.4 mag in the F330W filter. By contrast, the corresponding 
brightnesses decreased only half as much in 2003. It is clear from the timing of the observations 
that this difference is not simply due to a missed minimum in 2003.5. In order to have a 1-mag 
amplitude undetected in the 2003.5 observations, the 250 nm brightness would need either a strong 
negative "spike" only a few days long, or an appreciably delayed minimum in terms of phase in the 
5.5-year cycle, or both. Either of these possibilities would be as remarkable and unexpected as a 
major difference in amplitudes. (Unfortunately the minimum near 1998.0 is almost unknown.) 



Fernandez-Lajus et al.l (|2010l ) observed a similar but much less dramatic difference between 



events at visual wavelengths. In 2008-2009 they observed decreases of 0.15 to 0.25 mag in the 
UBVRI bands, only 0.02-0.03 mag deeper than during the 2003.5 event. It is unclear whether 
the difference in the depth of the minimum was small because they observed longer wavelengths, 
or because they measured the brightness of star plus inner ejecta rather than the star alone. 
Unfortunately we could not obtain suitable HST data at blue to red wavelengths in 2005-2009, 
because even the shortest integration time allowed for WFPC2 would have greatly over-exposed 
the star with any of the filters appropriate for photometry at A > 360 nm. 

The ground-based visual li ght curves suggest a minimum in 2009 between t = — 2 d and +3 
d (IFernandez-Lajus et al.1 120101 . observations at MJD 54858-54863; t was defined in Section 12.31 
above). In HST data, observations at t = —5 d and +12 d showed the lowest brightness, in 
agreement with their data. Other aspects of the 2009 event are consistent with broad-band optical 
and infrared observations of past events (jFernandez-Laius et al.l 120031 ; Ivan Gender en et al.l 120031 : 
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Whitelock et al.ll2004l ; iFernandez-Lajus et al.ll2010l ). The star began brightening at an accelerated 
rate about 40 to 50 days prior to the onset of the event. It brightened about 0.1 mag in F250W 
and 0.2 mag in F330W before it steeply declined. Then, after about 150 days in F250W and about 
100 days in F330W, the star returned to pre-event levels and resumed its long-term brightening 
trend. 

Why have the events become successively deeper in UV photometry (Figure [T])? Here is one 
line of reasoning. The primary star's ill-defined UV photosphere is located in the opaque wind, 
far outside the star's surf ace. If the wind den s ity is gradually decr easing on a time scale of ~ 
10 years (Section 1 above, iMehner et alj|2010bl . iMartin et alj|2006bl ). then the photosphere seen 
at most times must be shrinking and becoming hotter. During a spectroscopic event, however, 
the mass-loss rate may be enhanced, at least at low latitudes. This possibility depends on the 
secondary star's tidal / radiative influence and on the size of the primary wind's acceleration zone 
(|Davidsonlll997l . 1 19991 : Iffillier et allboOll : ISmith et allbood : IMartin et aDl2006ah . During an event 
the UV photosphere may therefore depend on different parameters than it does between events. 
(In principle, near periastron the secondary star might even be able to induce an outflow that is 
unrelated to the primary's usual wind. We do not advocate such an extr eme model here, but it 
illustrates the basic idea.) Meanwhile the Fe II forest around A ~ 260 nm ([Davidson et al.lll999bl : 
is sensitive to temperature and density. For these reasons, one should not be 



Hillier et al 



2001 



surprised if the UV photospheric size, temperature, and brightness during events follow a different 
long-term trend than they do at other phases of the cycle (Figure [1]). The same considerations 
apply to latitude an d longitude depend ences of the inner wind, even if there is no enhancement of 
total mass- loss rate (jSmith et al.ll2003l ). 



4. He II A4687 and X-rays 



He II A4687 provides extraordinary clues to the nature of the spectroscopic event. This is by far 
the highest-excitation feature known in rj Car's UV-to-IR wind spectrum, it ap pears only briefly at 
a cert ain stage in each event, and it p robably signals a floo d of very soft X-rays. ISteiner Damineli 
(|2004l ) first drew attenti on to it, but IMartin et al.l (|2006al ) disagreed with their flux measurements 
and interpretation, and iTeodoro et al.l (|201ll ) concur with the latter in most respects. Based on 
the Martin et al. analysis, we can summarize the relevant physics: 



• Since the observed feature almost certainly results from He ++ — > He + recombination, a tem- 
porary source of He + -ionizing photons (hv > 54 eV) is required. Shocked gas flowing through 
the wind-wind collision zone does not produce enough A4687 emission via normal cooling. 

• Nearly all authors agree that the two stars produce very little radiation above 54 eV. Therefore 
the relevant photons are most likely 54-500 eV X-rays produced in the wind-wind shock 
structure. Shocked gas of the primary wind, with pre-shock velocities below 600 km s -1 , is 
favorable for creating soft X-rays. 
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• The most suitable locale for A4687 emission is in the primary wind just before it encounters 
the colliding-wind shocks, and/or in locall y cooled condensations within the shoc ked region. 
Since the primary- wind shock is unstable (jPittard Corcoranl 12002 : ISokerl 120031 ) , these two 
choices may co-exist in roughly the same large-scale volume if one smooths over the complex 
small-scale structures. In either case A4687 is excited via photoionization by the soft X-rays 
mentioned above. 

• The most plausible energy source is the primary wind. A naive assessment predicts that the 
soft X-rays are inadequate for this purpose by a factor of 3-10, but several effects improve the 
efficiency. Martin et al. described radiative excitation effects that amplify A4687; instabilities 
in the shocked region tend to increase the number of very soft X-ray photons; and a brief rise 
in the primary wind outflow (hinted by other observations) would also help. With reasonable 
parameters, these details can enhance the He II A4687 flux by a sufficient factor. 

• The supply of soft X-rays can temporarily rise to very high levels if the fast secondary-wind 
shock becomes unstable like the primary-wind side. In that case the entire wind-wind interface 
"disintegrates" and "collapses," and a chaotic ensemble of subshocks and oblique shocks may 
exist for a few days or weeks. This phenomenon may explain the brevity of the A4687 flash 
as well as the disappearance of 2-10 keV X-rays. 



So me of these statements are controversial, but no other quantitative analysis has been pub- 
lished. ISteiner &: Daminelil (J2004J) proposed that He II emission occurs in the acceleration zone of 
the secondary wind, a much smaller region than those mentioned above. If one employs consistent 
physical parameters, their model predicts a A4687 flux two or thr ee orders of magnitude too weak 
(jMartin et al.1 l2006al : ISoker &: Beharl 12006 ) . ISoker Beharl (|2006l ) also focused on the inner wind 
of the secondary star, but their scenario was very different, emphasizing a collapse of the shock 
structure followed by accretion onto the secondary. They gave qualitative arguments for enhanced 
A4687 emission in specified circumstances, but did not quantify the excitation physics . Their mode l 
includes some appealing components which we note in later sections below. iTeodoro et al.l (|201ll ) 
recently indicated agreement with most of the above outline, and did not attempt a new theoretical 
investigation. In summary, the Martin et al. assessment is still the only detailed account of the 
radiative processes; and so far there have been no strong arguments against it apart from geomet- 
rical details. Below we shall mention various agreements between authors, and then some crucial 
disagreements. 
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4.1. The first He II Maximum during the 2009.0 Even 

Figure [2] shows a time sequence of the He II A4687 profile during the 2009.0 event, based on 
Gemini GMOS data. One can liken its appearance to a wave that first moves leftward in the figure, 
and then is reflected back toward the right. (Note the reversed asymmetry of the profile at t = +21 
d vs. —20 d. Here "reflection" alludes to the line profile, not a real physical reflection.) The greatest 
source of unc e rtainty in this feature's strength is the underlying continuum level, which, following 



Martin et al 



(2006 aj), we e stima te by interpolation between A4605 and A4744 Af£| Observations 



reported by lTeodoro et al.l (|201ll ) agree well with ours. At its maximum, the He II A4687 emission 
extended across 20 A or more (AV > 1200 km s 1 ) and had a much larger flux than one might 
guess from the apparent size of the line profile. This is true because the apparent "continuum" on 
each side of the line includes considerable A4687 emission; see Figure 3 in Martin et al., Figure 2 
in Teodoro et al., and the marks on the right side of our Figure [2J Part of the 20 A line width may 
conceivably be due to Thomson scattering rather than bulk velocities. Here we shall ignore some 
minor anomalies, e.g., the profile seen at t = —35 d was blueshifted less than either of its neighbors 
at -42 and -29 d. 

Figure shows the time development of equivalent width and Doppler velocity for He II A4687 
in our Gemini data. Here "equivalent width" refers to flux between 4675 and 4694 A (—770 to +450 
km s _1 ); other spectral features hide the farther line wings. Our equivalent width measurements 
are listed in Table "Velocity" in Figure [3] refers to the line's peak. As Teodoro et al. have 
noted, the main, double-peaked, 40-day A4687 flash in 2009 closely resembled the 2003.5 event. 
Its flux grew concurrently with the decline of 2-10 keV X-rays (Section 14.31 below); the maximum 
equivalent width agreed with the Martin et al. value for 2003.5 within measurement uncertainties; 
the time of significant brightness extended over 4-6 weeks; and the decline occurred in only one 
week. 

The most precise time marker for A4687 emission is the midpoint of its decline, which oc- 
curred at MJD 54843 or t = —17 d. The corresponding time in 2003 was close to MJD 52821 



(jSteiner &: Daminelill2004l ). and the difference of 2022 ± 2 d matches the consensus 2023-day spec- 
troscopic period (see Appendix). The decline-midpoint in our HST photometry at A ~ 250 nm 
occurred within a day or two of the same time (Section [3]) . Moreover, within measurement errors, 
the most negative A4687 Doppler velocity also coincided with the flux-decline-midpoint. (Steiner 
and Damineli noticed the same coincidence in 2003, and considered it evidence for an eclipse.) At 
the time of maximum A4687 brightness 8 days earlier, the line peak had Vdoppier ~ —310 km s _1 
(heliocentric); but then it rapidly moved to —420 km s~ 1 at the decline midpoint. Later we shall 



4 Throughout this section we use a time measurement t, with t — corresponding to MJD 54860, see Section \2. 31 
and Appendix. 



Steiner fc Daminelil l|2004f ) greatly unde restimated the peak A 4687 luminosity in 20 03 because they chose too 



high a continuum level, shown in Figure 3 in lMartin et al.1 (|2006af ). iTeodoro et al.l (120111 ) later adopted the Martin 
et al. method. 
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indicate reasons why the similarity between the 2003 and 2009 records is somewhat surprising. 



4.2. The "Second Episode" 

A second, smaller He II A4687 flare occurred several weeks after the first, around t ~ +20 d 
in Figures [2] and EJ Its significance will appear in Section 14.31 below; it did not match the 2003 
spectroscopic event. Independent of that question, two aspects of the "second He II episode" are 
noteworthy. First, the feature began to grow soon after its minimum at t ~ —10 d, and continued 
to do so for about a month. Therefore the emitt ing region was not entirely eclipsed by the primary 



wind as advocated bv lSteiner &: Daminelil (|2004j). A second fact is the rapidity of change in Doppler 
velocity. The first maximum or episode attained Vdoppier ~ —420 km s" 1 , but the emission that 
reappeared several weeks later had Vd ppi er ~ km s . Meanwhile the line profile's asymmetry 
reversed as mentioned earlier (Figure [2]). The overall velocity range greatly exceeds the maximum 
projected orbital velocity variation for either star with any proposed 5.5-year orbit. On the other 
hand, the line-of-sight wind velocities and post-shock velocities can easily differ by more than 400 
km s _1 at locations near the path followed by the secondary star near periastron. This statement 
may have interesting and controversial implications for the orientation of the orbit, see Section [4.41 
and Section [H 

Unlike the first He II maximum, the second episode in 2009 apparently differed from the 2003.5 
event. Here we disagree with Teodoro et al. about the 1992-2003 record, not 2009. The most 
reliable dat um for this question was the HST STIS observation at MJD 52852 in 2003, carefully 



assessed by iMartin et al.1 (l2006al ). Its timing corresponds to t = +15 d, only a few days before 
the phase of the second-episode peak seen in 2009. The 2003 STIS data showed no detectable line 
profile at that time, see Figure 5 in Martin et al. If one defines the emission feature as a bump 
or profile in the spectral tracing, then Martin et al. found a 3a upper limit of only 0.15 A for 
the equivalent width at MJD 52852. The points plotted in Figures [3H5] represent a different type 
of measurement, integrating a flux difference /(A) — / C (A) over a 19 A interval. Here / C (A), the 
underlying continuum level, is the main source of error. The 2003 STIS data point at t = +15 d in 
Figured] may therefore be an overestimate, but the absence of a bump in the spectral tracing (see 
above) indicates that it is not likely to be an underestimate. The observation is consistent with 
EW ~ zero at that time, and excludes EW > 0.30 A; see Sections 5 and 6.3 in Martin et al. If 
the 2003 event had a second He II episode matching the 2009 case, then the equivalent width at 
MJD 52852 should have been close to 0.9 A, which would have been obvious in the STIS spectrum. 
Conclusion: The 2003 event did not include a second He II episode matching that seen in 2009. If 
one did occur at the same phase, then in 2003 it must have been much weaker then 2009, even 
though the first He II maximum was practically alike in both events. But there is no good evidence 
for such an occurrence, and in Section [4.31 we favor the simple view that either 2003 had no second 
episode or else it was delayed. 



The ground-based data are extremely unclear on this question, especially since 



Teodoro et al 
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(|201ll ) disagree with ISteiner fc Daminelil (|2004l ) on major results even though they used similar, 
overlapping data sets. Only one data point is strongly pertinent: Teodoro et al. estimate EW(A4687) 
« 0.65 ± 0.16 A at MJD 52841, eleven days before the 2003 STIS observation noted above. In 
order to accept both this measurement and the STIS result, one would need to postulate a sudden 
decrease between them. This may be possible, but in that case the record would not resemble the 
2009 second episode. A simpler explanation is that the MJD 52841 datum was an overestimate (see 
below). Teodoro et al. also note a data point with EW w 0.56 A at MJD 48825 in 1992 (t sa +34 
d), and perhaps two weaker ones, but the doubts expressed below apply even more forcefully to 
them. 

Here we must digress on comparisons between ground-based data and HST. All ground-based 
spectrosco py of r? Car is seriou sly (not just mildly) contaminated by emission-line ejecta at r ~ 
Cf.'2 to 2" (jDavidson et al.lll995l ). The size of this effect depends on the spectrograph aperture size 
and on variable atm ospheric seeing. Mo reover, it was worse in the 1990's because the star was 
relatively faint then (IMartin et al.ll2006bl ). There is concrete evidence that narrow emission lines 
in the ejecta tend to mimic a weak "He II A4687" feature when the data are smoothed. Details can 
be found in iMehnerl (|201ll ); here is a brief summary. Gemini data and slightly decentered STIS 
spectra show many weak narrow spikes near 4687 A. Most of them are scarcely distinguishable 
from noise, except that their wavelengths recur in independent observations. This fa ct is not 
surpr ising, since the ejecta are known to produce thousands of weak unidentified lines (jZethson 



2001 



). When the spectra are heavily smoothed, some of these tiny features blend to produce a 
"A4687" bump with height ~ 0.5-1% of the continuum level. This bum p consistently appea r s nea r 
the same waveleng t h, an d it closely resembles the feature reported by ISteiner &: Daminelil ((2004) 
and lTeodoro et ah ( 2011 ) in data from Pico dos Dias (OPD) at times when A4687 was weak. The 
OPD feature tends to be stronger by a factor of ~ 2, which can be ascribed to worse atmospheric 
seeing, a larger spectrograph aperture, and th e fact that the apparent brightness ratio (e jecta/star) 
was much larger in the 1990's than it is today ([Martin et al.ll2006bl : iDavidson et al.lll995l ). (Despite 
this long-term trend, however, Teodoro et al. found the feature in data as late as 2010.) In summary, 
the OPD data show symptoms of extraneous emission- line contamination, at a level that is not 
surprising but which strongly affects most of the A4687 observations. HST STIS, by contrast, had 
good enough spatial resolution to separate most of the ejecta. We conclude that pre-2004 ground- 
based data have very low weight compared to STIS in 2003 and many ground-based observations in 
2009. Consequently there is no significant evidence for a second He II A4687 flare around t ~ +20 
d in 1992, 1998, and 2003. 

There is an interesting puzzle in the obser yational record. If w e assum e similarity among the 
1992, 1998, and 2003 spectroscopic events, then lSteiner fc Daminelil (|2004l ) depicted a conspicuous 
second He II episode at a significantly delayed time t. Their two relevant OPD observations occurred 
at t ~ +30 d in 1998 and +50 d in 1992, plotted in Figure [4] as green triangles. These appeared 
to be stronger than the most doubtful cases mentioned above. If valid, they imply second He II 
flares in 1992 and 1998, but at phases roughly 30 days later than the 2009 case - which would make 
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se nse for reasons n oted in Section [4.31 below. However, they contradict some 1992 ESO data listed 



by iTeodoro et al.l (|201ll ). and they are not mentioned in the latter paper even though it includes 
other data from the same instrument. Either there were delayed second episodes in 1992 and/or 
1998 (Steiner & Damineli) or not (Teodoro et al.). These inconsistencies reinforce our opinion that 
ground-based spectroscopy of rj Car in those years was far less robust than the 2003 HST results 
and various observations in 2009. 

In summary: One of the high-quality HST STIS observations in 2003 showed that the event 
in that year did not include a second He II episode matching the 2009 record; and this assertion 
is not contradicted by any reliable ground-based data. A weaker and/or later second episode may 
have occurred in 2003, as we discuss in Section [4.31 



4.3. He II, X-rays, and Shock Breakup 



Martin et al .1 (|2006al ) emphasized that He II A4687 was anti-correlated with rj Car's observable 
2-10 keV X-rays during the 2003.5 event; the A4687 flash occurred as the X-rays declined (Figure 
I3J). This fact appeared very consistent with a shock breakup as outlined below, and the 2009 
"second episode" now strengthens the case. 



Davidson! (12002 ) remarked that known shock instabilities, rather than the eclipse scenario 



that was popular at that time, can best explain the rapid disappearance of r\ C ar's 2-10 keV 



X-rays during a spectroscopic event. ISokerl (j2003l ) noted quantitative details, and Martin et al 



(|2006al ) emphasized the relevance of He II A4687 to this con cept (see Section 9.2 o f that paper). 



Various researchers later adopted essentially the same idea (IDamineli et al 



2009 



Teodoro et al. 



2011 



2008a 



Parkin et al 



Two specific variants, with different causes but similar consequences, 



have been proposed. 



1. A sho ck structure becomes unstable if radiative cooling exceeds expansion cooling (jStevens et al 



2. 



19921 ). The slow primary-wind shock of rj Car is very unstable in this re gard, but the faster 



secondary-wind s hock stabilizes the overall structure in calculated models (jPittard fe Corcoran 



20021 ; ISokerl 120031 ) . The secondary shock may become unstable near periastro n, causing the 



entire shock structure to disintegrate on a timescale of 10-30 days (Section 9.2 in lMartin et al 



2006a|) 



Soker Sz Behar ( 2006 ) drew attention to another phenomenon that Stevens Sz Kallman (jl99ol ) 



had studied for X-ray binaries in general. Near periastron, soft X-rays from the shocked region 
can alter the ionization state of the secondary star's wind. A higher degree of ionization tends 
to weaken the line-dr iven accelerat i on, re sulting in a slower wind speed ("radiative inhibition," 
a term later used bv lParkin et al.l (|2009l )). This in itself would reduce the 2-10 keV flux; but 
another consequence is that the balance of wind momenta is altered, pushing the shocks 
closer to the secondary star. In an extreme case the primary wind can entirely suppress the 
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secondary wind. 

Other instabilities certainly occur, such as Kelvin-Helmholtz which mixes gas from the two winds, 
and obvious thermal instability as shocked gas cools below 10 6 K; but the two processes listed above 
have mutually similar large-scale consequences and they might even develop together. Observable 
2-10 keV X-rays rapidly and tremorously disappear as the highest temperature decreases; a flood 
of soft X-rays is created by the chaotic ensemble of local shocks as the overall structure breaks up; 
and the secondary wind may temporarily cease to exist. (Alternatively, it might survive in a slower 
form). Martin et al. noted that empirically, without referring to physics, this view is supported by 
the growing unsteadiness of 2-10 keV X-rays around their maximum. 

Realistic simulation of these phenomena for rj Car would be extremely difficult, and many 
sets of parameters must be explored. Pending such models, one can estimate the relevant time 
scales from the overall size of the shocked region and local velocities within it: say 1-3 AU and 
100-800 km s _1 respectively. These values suggest At ~ 2 to 50 days, consistent with the observed 
behavior. 

Shock disintegration or collapse was originally only one among several competing explanations 
for the 2-10 keV X-ray crash during a spectroscopic event. Then, in 2003 the previously unrecog- 
nized He II A4687 feature behaved just as the shock breakup idea would have predicted. It appeared 
approximately when the hard X-ra y s peak ed, it grew as they declined, and then ceased abruptly 



after a few weeks. As iMartin et al.l (I2006al) noted, one wo uld not expect this sequence and timing 



in an eclipse model like that of lPittard Corcoranl (120021 ) . for example 



Teodoro et al.l (|201ll ) have recently advocated a shock breakup model that closely resembles the 
earlier conclusions of Martin et al. We differ from Teodoro et al. in two important respects: First, 
they state that He II A4687 originates "4-5 AU downstream from the apex" within the shocked 
gas. However, as implied in point 3 of the physics summary near the beginning of this section, 
there is no clear reason for this assumption. Wherever very soft X-rays occur, they indirectly create 
A4687 in the nearest suitable gas (which probably does not need to be shocked material). Most 
of the applicable energy supply for soft X-rays near periastron is found at distances considerably 
less than 4 AU from the pre-breakup apex or vertex. Moreover, it seems possible or even likely 
that soft X-rays and A4687 emission may continue for 5-20 days after the original large-scale shock 
structure has ceased to exist. Our second difference is that Teodoro et al. give little attention to 
the evidentiary role of the unforeseen He II second episode in 2009; see below. 

The 2009.0 event provides a valuable new clue, the He II second episode discussed in Section 
14.21 Figure H] shows that its behavior was a reversal of the first episode seen 40 days earlier. The 
A4687 emission reappeared and peaked around t ~ +20 d while the 2-10 keV X-rays were still weak; 
then it declined concurrently with the growth of the X-rays. The overall time scale was comparable 
to that seen in the first episode. In the shock-instability scenario, of course, we interpret this 
pattern as the re-formation of a large-scale shock structure when the relevant densities become 
sufficiently low for it to be quasi-stable. This view raises several questions. 
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As Figured] shows, the hard X-rays reappeared earlier than expected in 2009, around t ~ +30 
d rather than t ~ +60 d as seen in 1998 and 2003 (jKashi &: Sokerl l2009al ; ICorcoran et al 



20ipJ). Given this difference, should the earlier events have shown the same He II behavior 
pattern as 2009? The shock breakup hypothesis "predicts" He II A4687 flares around t ~ +50 
d in 1998 and 2003, i.e., just before the hard X-rays reappeared. Equally important, a 
A4687 outburst as early as the 2009 case, t ~ +20 d, would not be expected. The second 
He II episodes would very likely have been weaker in those earlier events, because the A4687 
production rate depends strongly on local densities. Consider, for example, a case with 
orbital eccentricity 0.9 and periastron at t = 0, fairly consistent with most proposed orbit 
models. Then the separation between stars was about 3.5 AU during the pre-periastron 
A4687 maximum (t ~ —20 d), but more than 6 AU at t « +50 d just before the X-rays 
reappeared in 1998 and 2003. If we assume an unchanged primary wind, these values imply a 
factor-of-three lower density at the expected time of reappearance. According to the Martin 
et al. analysis, this would reduce the A4687 flux by a factor of 3-10 compared to its pre- 
periastron maximum - depending on local gas geometry, velocity gradients, etc. Moreover, 
the primary- wind density may have been enhanced in the weeks before periastron (IDavidson 
19991 : iMartin et al.l l2006al ) . One should therefore expect the second He II episode to have 
been weak in 1998 and 2003, possibly undetectable. In 2009, on the other hand, the second 
episode occurred earlier with higher densities; which is the main reason why it is a valuable 
clue. (As we explain later, this statement includes a subtle complication involving a secular 
decrease in wind densities.) 

Is the observational record consistent with these expectations? The answer is "yes," but 
unfortunately the pre-2009 data were too sparse for confident conclusions. We noted in 
Section PO that the second expectation, no substantial A4687 flare at t ~ +20 d in 2003, is 
confirmed by STIS data. The question of a later flare , however, amounts to a c onflict between 



observ ations reported by ISteiner &; Daminelil (|2004l ) and those reviewed by iTeodoro et al 
(|201ll ). The former data strongly imply a A4687 episode peaking around t ~ +50 d in 1992 
and 1998, see Figure H] and Section [4.21 above. Teodoro et al., on the other hand, omit those 
two data points and show no definite A4687 emission after t ~ +35 d in 1992 and 1998, and no 
data points in that time interval in 2003. Neither of these alternatives would contradict the 
shock breakup idea (see above), but it would be useful to know whether or not a detectable 
second episode occurred. Since the OPD data are not public (unlike HST and GMOS), it is 
difficult to assess these alternatives. Meanwhile, the 2009 data and 2003 STIS data constitute 
the only satisfying results on this question. 

Why did the He II emission and hard X-ra ys reappear earlier than expected in 2009? Referring 
only to the X-rays, iKashi Sokerl ()2009al ) proposed an explanation based on a decrease in the 
primary star's wind. They had a specific model in mind, but the basic idea works for others 
as well. Since 1999 there have been hints that t? Car's wind density is (or was) becomin g 
less dense on a timescale of ~ 10 years (e.g., IMartin et al.1 l2006bl : IDavidson et al.l 120051 ). 
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Mehner et al.l (|2010bl ) recently reported strong spectroscopic evidence for such a change. 
Qualitatively, at least, one expects the wind-wind shock structure to recover sooner after 
periastron if the primary wind is less dense; because in that case the worst instabilities are 
weakened. At first sight the 2009 He II second flare may seem paradoxical in this view, 
because, as noted above, lower densities tend to reduce the A4687 emission efficiency A 
partial explanation is that the shock-region density at the time when it re-forms is not the 
same as the primary wind density measured at some fixed radius; the two stars were closer 
together at t ~ +20 d in 2009 than they were at t ~ +50 d in 2003. This brings us to the 
next question which also concerns timing and separation between the stars. 

Suppose, as stated above, that the primary star's wind was less dense in 2008-2009. If so, 
why didn't the shock structure survive later than expected before periastron? With reduced 
density p w i n d( r )> one would expect the shocks to disintegrate at a smaller star-star separation, 
i.e., later. In fact the 2009 X-ray crash and the main He II A4687 flash occurred almost 
exactly 2023 days after the 2003 event (Section 14.11 and Figure HJ. Therefore, perhaps the 
characteristic densities at that time did not directly represent the primary star's "normal" 
steady mass-loss rate. Instead, the inner wind may have been tidally or radiatively enhanced 
during the weeks before periastron, causing the relevant density to be approximately the same 
for each spectroscopic event. This hypothesis is not arbitrary or ad hoc, since it also helps to 
explain a photometric puzzle m entioned at the end of Section [3] above, and it was suggested 



long ago on other grounds (e.g.. iDavidson 



1999 



Martin et al.ll2006al ). 



Teodoro et al.l (1201 If) have express ed agreement with most, but not all, of the points that we 



have quoted from iMartin et al.1 (|2006al ) . One disagreement pertains to the He II "second episode" 
seen in 2009. Above we emphasized the anti- correlation between A4687 flux and hard X-rays at 
critical times. Teodoro et al., however, state instead that the two are correlated but A4687 is 
delayed by 16.5 days. We strongly disagree, based on semi-theoretical reasoning plus an empirical 
fact. (1) Those authors interpret the 16-day delay as the time required for shocked gas to flow from 
a favorable X-ray region near the vertex or apex of the umbrella-shaped shock structure, to another 
location where it has cooled enough to produce He II emission. For this explanation one needs a 
localized density enhancement, which must pass near the shock vertex and later reaches some 
particular radius at a well-defined time after flowing outward within the shocked zones. But these 
are unlikely assumptions in the context of an extended, chaotic, unstable shock structure. There 
is no clear reason for the localized density enhancement; only a tiny fraction of material flowing 
outward from either star passes close to the shock vertex; gas is vigorously mixed and diffused 
within the shock structure; the overall geometry is probably unstable at the critical time as noted 
above; and there is no reason to assume that He II emission occurs only beyond a certain distance 
from the vertex. In other words the proposed rationale is very unclear. (2) Our second, more 
empirical reason is simpler: The 16-day delay proposed by Teodoro et al. is obviously inconsistent 
with the A4687 second episode observed in 2009 (Figure E|) . Maximum He II emission occurred 
after the hard X-ray peak in the first episode, but before the X-ray recovery in the second. 
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In summary: Strictly speaking we cannot prove the shock breakup scenario, but the combined 
A4687 and X-ray data are impressively consistent with it. In terms of logic, the role of He II 
A4687 has been as follows. First, this feature was not recognized when the breakup idea was first 
proposed; but then it turned out to fit into that scenario in a remarkably natural way. Later, when 
a post-periastron "second He II episode" appeared in 2009 as discussed above, it too matched quite 
naturally, and the observed secular decrease in wind density provides a good explanation for the 
differences between 2003 and 2009. In neither case was there any need to alter the concept to fit 
new observations. 



In terms of theoretical development, iDavids on ( 200 2]) and ISokerl (12003!) stressed that shock 
breakup is a reasonable concept for 77 Car; then IMartin et al.l (|2006al ) and ISoker &: Beharl ((2006) 
argued that it is a very likely one. Parkin et al. ( 20091 ) later cited other X-ray clues, notably that gas 
dynamic computations have not sustained the competing eclipse scenario. Given these facts, and 
lacking a viable alternative explanation, we conclude that 77 Car's colliding-wind shock structure 
does, indeed, disintegrate and collapse during a spectroscopic event. Two separate questions are 
(1) whether the primary star has a mass-loss outburst at about the same time, and (2) whether 
the secondary star accretes material then. See Section [6] below. 

Is He II A4687 emissi on d etectable in 77 Car on ly during spectroscopic events? According to 



Steiner fc Daminelil (12004T) andlTeodoro et al.l (|201ll ). this feature has been present at other times. 



On the other hand, IMartin et al.l (l2006al ) found no broad He II in non-event HST STIS spectra of 
the central star. STIS mapping data in 2009 June and December show no broad A4687 emission in 
extended regions around the star. Analyzing observations made with HST STIS, Gemini GMOS, 
and Irenee du Pont B & C, we find that a weak emission feature exists but it is only about 3 A 
wide, i.e., narrower than the He II emission seen when it is bright. As noted in Section 14.21 it is 
most likely a blend of very weak, narrow, unidentified lines and not broad He II emission (jMehner 
20 lib . 



4.4. He II A4687 in the Reflected Spectrum from the Pole (FOS4) 



Spectra reflected by dust in the Homunculus nebula give surprising new information about 
velocities in the 2009 event. The known geometry of the bipolar Homunculus allows us to correlate 
each position in the southeast lobe with stellar latitude, assuming that the polar axis of 77 Car is 
aligned with the Homunculus axis. FOS4 is a location near the center of the SE lo be (the one nearer 



Davidson et al. 


1995; 


Zethson et al. 


1999) 


(Davidson et al. 


2001; 


Smith 


2006 


). FOS4 



75° latitude ( Smith et al J 12003 



•ur direct view probably represents ~ 45° latitude 
also useful for another reason, namely that ground- 
based spectra there are less contaminated by nebular lines than direct observations of the star are. 
(The reasons for this fact are not entirely clear, but generally speaking our direct view of the star 
appears to have more circumstellar extinction than the average line of sight.) FOS4 is located 
about 3'/7 south and 2 / . / 5-3 / /5 east of the central star. 
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A spectrum reflected in the Homunculus has a light-travel-time delay At corresponding to 
the additional path l ength, and an ex tra Doppler shift AV due to the "moving mirror" effect 



(jMeaburn et al,lll987l ; IStahl et al.ll2005l ). With standard assumptions about the Homunculus, these 
are related by AV/c = At /(age), where "age" means elapsed time since the reflecting material was 
ejected in the 1840's. The value of AV depends on location in a s traightforward, observ able way 



that depends on the shape of the Homunculus lobe, see Figure 4 in lDavidson et al.1 (|2001l ). In fact 
we can define FOS4 as the location where AV = +10 km s _1 , which implies At ~ 20 d for the 



2003 and 2009 spectroscopic events. IStahl et al.l (|2005l ). however, reported a time-delay of only 10 
days for appearance and dis appearance of He II A4687 in 2003, comparing VLTUVES observations 
of FOS4 to observations bv lSteiner &; Daminelil (|2004l ) which were centered on the star. If correct, 
this would cast doubt on standard analyses of the reflection model, Homunculus expansion, etc. 
The good time coverage of our Gemini GMOS data during the 2009 event, both on the star and on 
FOS4, allows us to re-evaluate this result. Figure [5] shows the equivalent width and radial velocity 
measurements of He II A4687 in spectra of the star in direct view and reflected at FOS4. It also 
includes re-measured equivalent widths in UVES data during the 2003.5 event, shifted by 2023 
days. The GMOS data show that the time-delay in the observed equivalent widths as well as radial 
velocities, between the direct view on the star and the pole-on view at FOS4, is about 18 days 
which confirms the expected At ~ 20 d within the attainable accuracy. Also, re- measurement 



of the UVES data during the 2003.5 event according to the method employed by iMartin et al 



(|2006al ) show that the earlier data are in accordance with the GMOS 18 days value. We detected 
a positional gradient of At, which agreed with the simple model within the measurement errors. 
The geometry of the reflection process therefore appears satisfactory. 

Unexpectedly, we find that the behavioral pattern of He II \4687 emission is very similar when 
viewed from different directions, i.e. in direct view of the star and reflected at FOS4, when we 
take the time-delay into account. Values for the equivalent widths and radial velocities are slightly 
smaller at FOS4 than in the direct view. Figure [5] shows these results, with Doppler velocities 
corrected for the moving-mirror effect AV. He I A4714 seems to exhibit similar behavior, but 
cannot be measured well be cause its line profile changes shape in a more complex way than He II 
A4687 (Section [5] below, and lMehnerlhoilh . These results suggest that velocities of helium lines are 
not simply related to orbital motion of the secondary star, if we assume that the orbit inclination 
is i ~ 40-45° like the Homunculus midplane, see refs. cited above. In standard models the view 
from FOS4 should be almost perpendicular to the plane of the orbit, and therefore no large radial 
velocity variations should be observed there. Thus the data at FOS4 are very surprising. 



As examples, let us mention just two of the published models. ISoker &; Beharl (|2006l ) proposed 
that helium lines originate in the acceleration zone of the secondary star. Assuming we view the 
orbit at an inclination of 40-45°, then the He II Doppler velocity variation observed in our direct 
view, ~ 400 km s _1 , implies values of ~ 550 km s _1 in the plane of the orbit. The variation at 
FOS4 is ~ 250 km s" 1 ; if we interpret this too as the projection of a velocity in the orbital plane 
(70-80° inclination at FOS4), then true values of ~ 750-1500 km s _1 are needed. In order to attain 
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high enough projected flow velocities, the Soker and Behar model would therefore require FOS4 to 
"see" He II A4687 emission formed farther out in the acceleration zone than the emission in our 
direct view. This is not easy to arrange, because the emission process is inherently isotropic. Even 
more significant, one would also need to explain velocity variations that greatly exceed any credible 
orbital velocity. We shall comment later on possible alterations of the assumed inclination. Similar 
arguments apply to models in which the helium emission originates in or near the wind- wind collision 



region; for instance the He I interpretation by iNielsen et al.1 (|2007l ) is quite inconsistent with the 



FOS4 data. No matter whether one ascribes the observed Doppler variations to orbital velocities 
or to successive illumination of wind regions near the moving secondary star, FOS4 was expected 
to differ from a direct view of the star. The differences should have been far more conspicuous than 
our observations show. 

Given the consistent time-delay at FOS4, it is very difficult to imagine that the above result is 
mistaken; but here we can offer only conjectures to explain it. One might reconcile the data with 
a different inclination i, such that the projected orbit appears alike from both our point of view 
and FOS4. But that idea would "multiply the hypotheses" by requiring two unlikely assumptions: 
(1) The orbit plane must be tilted 20-25° from the Homunculus midplane, and (2) the azimuthal 
direction of the tilt must be aligned with both our line of sight and that of FOS4, so they share 
nearly the same projected velocities in phase with each other. In other words, both the latitude 
and longitude of the tilt must be suitable. This recourse seems too artificial and ad hoc to be 
an appealing first choice. A better explanation might be that the observed Doppler variations 
represent "global" changes in outward velocities, roughly spherical outflows that are not given a 
strong directionality by the secondary star. (An extreme form of this view might even return to the 



single-star mass ejection proposed by lZanella et al.lll984l . We do not intend to go that far, but the 



question is worth contemplating.) A shock breakup model may conceivably act in a quasi-spherical 
way, with chaotic random velocity components during the critical time. The same statement applies 
to enhanced mass-loss from the primary star. Independent of these thoughts, our FOS4 results cast 
doubt on attempts to derive orbit parameters from apparent emission-line velocities. 



5. The Changing Wind Structure during the 2009 Event 

HST STIS data revealed that some spectral features depend on viewin g direction and that 



the global stellar wind geometry changes during the cycle (jSmith et al.l l2003). The most dramatic 
effects occur at low latitudes, while the dense polar wind remains relatively undisturbed during an 
event. Departures from spherical symmetry are critical for theories of winds and instabilities in the 
most massive stars and we therefore re-examine selected spectral features at differing latitudes in 
our Gemini GMOS data. Smith et al. analyzed only three epochs after the 1998 event; 1998 March 
(phase = 0.04), 1999 February (phase = 0.21), and 2000 March (phase = 0.40). The improved 
time-sampling of the GMOS observations makes it possible to monitor changes before, during, and 
after the 2009 event. Our results do not all agree with Smith et al., and some of the differences 
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may signal real changes. 



Hydrogen: Ha and H/3 emission lines are so bright in rj Car that all Ha and many of the H/3 
observations centered on the star were overexposed in Gemini GMQS obse rvations 1^1 For this objec t 
H7 is usually contaminated by other emission lines, so we analyzed H<5. iRichardson et al.l (|2010i ) 
have described the behavior of Ha in 2008-2009, based on a large number of observations with 
higher spectral resolution; but Ha samples differ ent properties of the system because it originates 
at much larger radii than H.8 ^Hillier et alJboOlh . 



Figure [7] exhibits variations observed in the H<5 P Cyg profile at several latitudes. This figure 
shows spectral tracings seen at the star and four locations to the southeast, at several phases close 
to the 2009 event. GMOS observations before the event, from 2007 June to the beginning of 2009 
January, revealed no significant changes in the H<5 P Cyg profiles. At those times (t = —353 and 
—82 d), substantial P Cyg absorption was observed only at higher latitudes. This fact is usually 
considered evidence that the de nsity and /or ioniz ation structure of ry Car's current stellar wind 
outside an event is nonspherical ([Smith et al.ll2003l ). For most of rj Car's spectroscopic cycle, wind 
densities are expected to be higher n ear the poles, in accordance with theories of equatoria l gravity 
darke ning i n massive rotatin g stars (jMaeder &: Mevnetl 120001 ; iMaeder &, Desiacqued l200ll ; lOwocki 
20051 ). But IRichardson et al.l (|2010l ) have noted an alternative explanation, in which the secondary 
star prevents hydrogen P Cyg absorption in the primary wind. A basic obstacle to settling this 
question is that the fractional population of H° in the n = 2 level, essential for Balmer absorption, 
is an intricate theoretical question that has not yet been explored for 77 Car's wind. 



Standard theory predicts higher wind velocities at high latitudes, and ISmith et al.1 (120031 ) 
reported such a result extending almost to ~ —1000 km s -1 in the 2000 March STIS data (see their 
Figure 5). However, we find Voo ~ 500-550 km s _1 for H<5 absorption at all latitudes during 77 Car's 
normal state (Figure [7]). One reason for this discrepancy may be that we use a different method to 
align the spectra. Smith et al. corrected for AV in the expanding nebula (Section 14.41 above) by 
aligning the blue side of Ha; but that reference point in the line profile is itself affected by details 
of the P Cyg absorption. We use, instead, several forbidden lines that are known to originate in the 
Weigelt knots with constant velocities much smaller than the discrepancy in question. Also, the 
velocity structure of the wind may change throughout the cycle; Smith et al. used observations at 
phase 0.40 while our relevant observations are at phase 1.83 and 1.96. Our GMOS observation at 
phase 2.17 cannot be used to investigate this issue, since the H5 profile had obviously not returned 
to its normal state. H5 should be practically as good as Ha and H/3 for this purpose, since its 
absorption component is stronger, relative to emission, than in Ha. A real change in the velocities 
and latitude structure may have occurred between 2000 and 2008. 



As already discussed by Smith et al. ( 2003 ) and Stahl et al. (j2005l ). the strong latitude depen 



6 Gemini GMOS does not support exposure times below 1 s. 
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dence of Balmer P Cyg profiles does not apply during the events. In only a few days, between 
t = —25 d and —20 d, the P Cyg absorption at lower latitudes appeared and strengthened to the 
same depth as at higher latitudes. Strong absorption at all latitudes was observed until t m +47 d, 
i.e. for at least 70 days. Observations from t = +89 d to +344 d showed only weak absorption at 
low latitudes while high latitudes continued to have strong P Cyg profiles; the system had almost 
returned to its pre-event state. 

P Cyg absorption is always present in the higher Ba lmer lines at all latitudes, consistent with 
their formation regions closer to the star (|Weis et al.l l2005). Still, during the events, their absorption 
deepened at lower but not higher latitudes. 

Doppler velocities of Balmer l ines are difficu l t to as sess because they may include two different 
but unresolved parts. As noted in lMehner et al.1 ( 2011 ). the main velocity of H<5 emission remains 
fairly steady, while a second component appears to vary like He I (see below). The two components 
overlap so much that neither can be studied individually. 



He I: Both the emission and the absorption components of He I vary in strength and radial 
velocity throughout the cycle, see Figure [6] for measurements on the star. The equivalent width 
of the He I emission was mostly constant during the cycle, increased before the 2009 event, and 
then dropped into a temporary minimum. However, the emission line never disappeared entirely as 
would have been expected if the event were a true eclipse. The emission lines shifted monotonically 
blueward throughout the cycle, terminating with an abrupt, large velocity shift of over —100 km 
s _1 to velocities of about —250 km s" 1 near the event, followed by a sharp rise to almost zero 
radial velocity - very similar to the He II A4687 behavior noted in Section 0] above. The radial 
velocity of the absor ption lines showed a similar pattern with velocity shifts between —300 km s _1 
and —600 km s ISmith et al.l ([20031 ) found that while the He I emission faded at low latitudes 
during the 1998 event, the emission was relatively undisturbed at higher latitudes. However, GMOS 
observations show that the equivalent width of the He I emission at FOS4 during the 2009 event 
followed basically the same patterns as the emission directly on the star. 

Figure [8] shows tracings of He I A4714 in GMOS observations at several latitudes close to 
the 2009 event. Outside the event, from 2007 June to 2008 July, He I lines had strong P Cyg 
absorption in spec tra at low latitude s and only very weak P Cyg absorption in spectra at higher 
latitudes (see also ISmith et al.l 12003). Because of their limited time-sampling, Smith et al. were 
not able to observe changes in the He I P Cyg absorption. This led them to conclude that for most 
of the cycle He I absorption is present on the star but not at higher latitudes. Their Figure 18 
implies that the absorption at low latitudes would disappear during the event, though this is not 
explicitly stated in the paper. However, in our data we find that shortly before the 2009 event, 
the He I absorption increased at higher latitudes. At t = —82 d the P Cyg absorption was strong 



7 lNielsen et al.l ( 2007T ) found similar values during the 2003.5 event. 
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at all latitudes. Then, over the next 2 months the absorption weakened at all latitudes. Between 
t ~ —5 d and t ~ +89 d almost no absorption was observed. By t ~ +176 d the He I profile 
had returned to its normal state displaying strong absorption at lower and almost no absorption 
at higher latitudes. Note that changes in the He I P Cyg absorption occurred at least 2 months 
earlier than those in the HIP Cyg absorption, and the return to the pre-event state took up to 3 
months longer for He I. The presence of He I P Cyg absorption at all latitudes s hortly before the 
event may have implications for the shell ejection model favored by lSmith et al.1 (J2003J), since this 
observation is not directly accounted for by that model. 

The geometrical volume filled by ionized helium around ij Car is highly relevant, since He I 
absorption as well as emission lines depend on He + — > He recombi nation. The zones of interest are 



photoionized chiefly by the hot secondary star, see iHumphreys et al.1 (|2008l ) . iMehner et al.l (l2010al ). 
and refs. therein. If t he primary wind between events has become less dense in the past few years 
(jMehner et al.ll2010bl ). then the He + zone photoionized by the secondary star should now occupy a 
much larger volume than it did at earlier times. A crude assessment of the helium-ionizing photon 
supply suggests that the He + zone most likely now wraps around the inner primary wind at times 
other than spectroscopic events. In other words, the pseudo-hyperboloidal He + /He° ionization 
front in the primary wind may now be concave toward the primary star. If so, then most lines of 
sight to the primary star - including that for FOS4 - must pass through some He + in the primary 
wind, unlike the case 10 years ago. Detailed work will be presented in a future paper. 



Fe II: Narrow lines and broad stellar wind emission of Fe II are strong in dir ect spectra o f p Car , 
but much fainter in the reflected spectrum at FOS4. As already noted by ISmith et al.l (|2003l ). 
the Fe II lines resemble Balmer lines in that the Fe II P Cyg absorption increases with increasing 
latitude and that the emission is weaker at higher latitudes. Figure [9] shows broad stellar wind 
emission of Fe II in spectra at FOS4. Spectra before t ~ —82 d showed only very weak absorption at 
FOS4, with maximum strength at ~ — 400-450 km s _1 , but the absorption feature then deepened 
and stayed strong until t rj +176 d, i.e. for about 250 days. The absorption was strongest around 
t ~ +10 d. Th e deep ening of Fe II absorption at FOS4 was also observed during the 2003.5 event 
by IStahl et al.l ([2005) . Other species, such as Cr II, Mg I, and Ti II, also developed absorption 
lines. 



Smith et al.l (J2003J) argued that Fe II lines are formed in the same outer wind regions as H I 



because these two species show similar latitude dependence and are likely to be coupled by charge 
exchange. However, their similarity is only true outside an event. During or close to an event 
the lines behave very differently; in contrast to H I, Fe II does not develop absorption in direct 
view. (H I lines develop strong P Cyg absorption at lower latitudes while the pole remains almost 
unchanged.) Theoretically, H I Balmer absorption is not really like most Fe II absorption, since the 
hydrogen n = 2 level has a much higher energy than most Fe II lower levels (10 eV compared to 
0-3 eV). The observed Balmer absorption lines may involve H + — > H° recombination, analogous to 
He I noted above. 
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6. Critical Differences in the 2009 Event — A Summary, Conclusions, and Problems 

6.1. Primary Results 

We monitored rj Car with HST WFPC2 and Gemini GMOS throughout its 2009 spectroscopic 
event. Good time coverage with the GMOS slit oriented at a constant position angle made it possible 
to monitor spectroscopic changes as seen from a range of directions, some of them via reflected light. 
In this paper we have described several important differences compared to previous events, some 
of them quite unexpected. These results lend strong support to one important concept, the idea 
of shock breakup and collapse near periastron. Regarding other aspects of the problem, the new 
information helps to exclude some proposed models, while - as usual for this topic - it also deepens 
some of the puzzles. 

HST WFPC2 observations show that the minimum in the UV was much deeper for the 2009 
event than for the 2003.5 event (Section [3|). One possible explanation involves a mass-loss outburst 
or similar phenomenon, discussed in Section [6.21 below. 

Contrary to expectations, we find that the behavior of the He II equivalent width and radial 
velocity reflected at polar location FOS4 are very similar to direct observations of the star. For 
this purpose an observed light-travel-time delay At f» 18 d and moving-mirror redshift AV ~ +100 
km s _1 must be taken into account. Since the observed At agrees well with the predicted value, it 
confirms that FOS4 really does "see" the star from a polar direction. The observed radial velocity 
behavior at FOS4 is a surprise, because, contrary to proposed models, it fails to show major 
differences compared to our direct view of the star. This statement includes other effects as well as 
the orbital velocity. 

We found that the "second He II A4687 episode" in 2009 was strongly anti-correlated with 
the X-rays, like a time-reversal of the main A4687 outburst seen in 2003 and 2009. In Section 14.31 
we argued that this result strongly supports the shock structure breakup/disintegration/collapse 
hypothesis, first proposed a decade ago as an explanation for X-ray behavior during a spectroscopic 
event. Strictly speaking this idea has not been proven, but now the following facts together have 
made it the most probable model. 



1. Likely physical conditions near periastron appe ar suitable for the relevant instabilities (jSokei 



20031 : Martin et al.ll2006al : ISoker fc Beharil200(Sl ) 



The He II A4687 outburst was not yet recognized when shock breakup was first proposed in 
2001-2003, see Section f4. 31 above. But after it was discovered and measured, this phenomenon 
turned out to be very well suited to the idea. 

The 2009 event featured a new development, the second A4687 episode about 30 days after 
the first. As we discussed in Section \4. 31 this was beautifully accordant with the earlier-than- 
expected reappearance of 2-10 keV X-rays. 
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4. The time-scales are reasonable as noted earlier. 

5. Prom a theoretical viewpoint, other proposed explanations (especially eclipses, which most 



authors favored un til two or three years ago) have much worse difficulties ([Martin et al. 



Parkin et alJl2009h . 



2006a 



But shock breakup is still a hypothesis, albeit a very strong one, not a quantitative theory. Soft 
X-ray production during the collapse, and structural recovery later, have scarcely been investigated 
yet. A broad range of parameter space must be explored. 

Let us emphasize that shock breakup would explain only some aspects of a spectroscopic event. 
Other phenomena, more fundamental for the basic physics of rj Car, are probably required in order 
to explain the spectroscopic and photometric changes. Some of them are noted in Section HT21 below. 

Results on spectral features of hydrogen, He I, and other species (Section [5]) are too diverse 
to summarize briefly. H ydrogen P C y g pro files at different latitudes throughout the cycle behave 



as already discussed in ISmith et al.l (|2003l ). Outside the 2009 event, HIP Cyg absorption is 
observed at higher latitudes but not at lower, while during the event HIP Cyg absorption is also 
strong at lower latitudes. GMOS data show that the P Cyg absorption at lower latitudes appeared 
suddenly within only a few days and was prese nt for at lea s t 70 days. We do not find a higher 



terminal velocity at higher latitudes as found by lSmith et al.1 ([2003). The exact reasons have to be 
addressed in the future. 

Helium P Cyg profiles showed an additional incident not discussed by Smith et al. Those 
authors found that outside the events He I P Cyg absorption is present at low latitudes but absent 
at higher latitudes, while during the events the absorption disappears at low latitudes, too. Our 
GMOS data show that shortly before the 2009 event He I absorption increased at higher latitudes to 
similar strength as at low latitudes. Then the absorption decreased slowly at all latitudes. Changes 
in the He I lines were observed already 2 months before the H I lines showed any changes and they 
returned to their normal "no-event" state up to 3 months later than the H I lines. 

Fe II absorption, only present at higher latitudes, became very strong during the 2009 event 
for several months and behaved different than the H I lines. 

Further analysis is required with regard to the cause of those observed latitude dependent 
changes throughout the events. Is a minor shell ejection sufficient to explain them? Or do changing 
ionization fronts in the primary wind caused by influences of the secondary star and/or the moving 
wind-wind shock structure play a role? 



6.2. Unsolved Problems 



A number of essential questions have not yet been answered even after years of observation 
and discussion. Each requires theoretical work that no one has attempted at a realistic level of 
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detail, and our 2007-2010 results are pertinent to some of them. 



Does the secondary star accrete material during a spect roscopic event? This pqssibil 



ity has been emphasized especially by Soker and colleagues (e.g. ISokerl 120031 ; ISoker Beharll2006l ; 
Sokerl 120071 ; iKashi fc Sokerl l2009al ) . As we noted in Section 14.31 near periastron the primary wind 
may entirely suppress the secondary wind, thereby allowing accretion onto the secondary star. This 
possibility can be separated from other aspects of those authors' model. It is a very appe aling idea 
because it may explain a long -standing paradox concerning r\ Car's He I emission lines ((Davidson 
19991 ; Humphreys et al.l 12008). These features depend chiefly on photoionization by the hot sec- 
ondary star, they weaken during a spectroscopic event, and they were not present before 19^1. 
Merely saying that "the secondary star was engulfed in dense gas" does not constitute an expla- 
nation, since helium-ionizing photons {hu > 24.6 eV) inevitably generate He I recombination lines 
in the primary wind even at densities far above normal. R oughly spea king, the brightness of He I 



emission measures 



the far-UV photon supplyjfl However, as Soker ( 2007 ) noted, accretion can lower 



the secondary star's effective temperature by slightly expanding its photosphere. Even a modest 



reduc tion in temperature substantially cuts the output of helium-ionizing photons (jMehner et al 



2010al and refs. therein). Humphreys et al. remarked that rj Car's mass-loss rate may have been far 
above 10 -3 M Q yr 1 in the years 1900-1940; the primary star's wind density near the secondary 
star then may usually have been as l arge as it is today nea r periastron - i.e., very likely enough to 
crush the secondary wind. (See also lKashi &; Sokerl l2009bl . concerning r/ Car's Great Eruption.) 



Accretion thus provides a satisfying explanation to the He I puzzle; but is it correct? Helium 
emission lines did not disappear in 2009 (Figure E|) , but on the other hand this was presumably the 
least dense spectroscopic event on record. Truly realistic simulations of three-dimensional accretion 
and resulting photospheric temperature will be extremely difficult, but they may prove necessary. 
At present we cannot be certain that the secondary wind is suppressed near periastron. 



Does each spectroscopic event include a mass-loss outburst b y the primary star? This 
suggestion arose lon g ago because eclipse models appeared inadequate ( Davidson! 1999 ) . and because 
Zanella et al.l (| 19841 ) had discussed the same idea in a single-star context. Today the shock breakup 
concept probably explains the X-ray and He II A4687 behavior, but other spectral features involve 
larger spatial regions in the unshocked primary wind. Throughout this paper we have mentioned 
indications that the primary wind may have been disturbed. Most of them have been noted before, 
and none is entirely satisfying, but their combination is highly suggestive: 



Our UV photometry showed much deeper minima in 2009 than in 2003. In Section [3] we noted 



Mehner et all (|2010aT ) estimated T off « 40, 000 K and 10 5 L < L < 10 6 L Q for rj Car's secondary star. We 



estimate that such a star does indeed produce enough helium-ionizing photon s to account for the observed equivalent 
widths of rj Car's He I lines. In addition to the reasoning use d for nebulae ( Ostcrbro ck fc Ferla nd 2006), one must 
include a special enhancement factor explained in Section 6 of Humphreys et al. ( 20081 ). 
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why this may be evidence for extra material in the system during the early stages of the event. 
Qualitatively, a higher density due to temporary causes would affect the Fe II forest and UV 
photometry in the desired direction. Shock breakup, by contrast, does not provide such 
an obvious explanation. A similar argument can be made regarding near-infrared free-free 
emission. 

Temporarily enhanced gas densities would help to destabilize th e shock structure, while also 



supplying extra energy for He II A4687, as iMartin et al.1 (|2006al ) explained. 



In Section 14.31 we noted that the 2008-2009 X-ray crash "should have" occurred later than 
predicted, since it is now known that the primary wind had gre atly decreased since 2003 
(jMehner et al.l l2010bl : iKashi fe Sokerl l2009al : ICorcoran et al.ll2010l ). Extra ejected material 
dependent on other factors, however, would have made the secular density trend temporarily 
irrelevant. 



In Section 14.41 we stated that a reflected polar view of the 2009 event closely resembled the 
behavior seen directly, including the Doppler variations. Any effect involving orbital motion 
should have appeared different between the viewpoints. A conceivable explanation is that 
velocities during the event may have represented instead a varying quasi-spherical outward 
flow from the primary star. 

Complicated changes in P Cyg absorption features, described in Section [5] above, suggest that 



colum n densities increase especially at low latitudes during an event (compare ISmith et al 
20031 ). 



If the Eddington limit is taken into account (and, perhaps, rotation), tidal force s are very 
likely strong enough near periastron to alter the primary wind-acceleration zone ([Davidson 
19971 ). Various instabilities may have broadened the flow. 



Of course we do not claim that a mass ejection event would explain everything. But it would help 
with all the above points, and there is no clear evidence against it. 

The outburst conjecture has both milder and stronger variants. For instance, the total mass- 
loss rate might remain constant during the event, but its la titude dependence briefly changes, 
causing densities to increase at low latitudes ([Smith et al.ll2003l ). At the other extreme, the base of 
the wind might be affected, not just its acceleration zone. That case would signal an undiagnosed 
surface instability in the star itself. How can this question be investigated? As a beginning, 
one might explore theoretically whether the photometric and spectroscopic observations can be 
explained by shock breakup without a primary wind outburst or disturbance. 

The primary wind de nsity has been decrea sing and may continue to do so until the star has 
a normal line-driven wind (jMehner et al.ll2010bl ). If so, then the wind-wind shocks will very likely 
become stable even at periastron, ending the series of X-ray and He II "events" as we know them 
today. This might even occur within the next two or three 5.5-year cycles. If, on the other 
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hand, the secondary star is capable of triggering a tidal/radiative outburst near periastron (as 
discussed above), then the X-ray events will continue to occur. In that case the X-ray luminosity 
will be diminished, but not as severely as the wind density. (This statement is based on simple 
considerations of the wind-wind momentum balance and cooling rate.) 



How can the orbital parameters be estimated? The FOS4 results bolster our opinion that 
observed Doppler velocities cannot be used for this purpose until they are much better understood. 
There have been attempts to derive orbit parameters by directly identifying r\ Car's observed 



Doppler variations with orbital velocities, like a classical spectroscopic binary (e.g.. iDamineli et al 



1997 



Nielsen et al.l 120071 ; iKashi Sokerll2008l ). This approach is questionable because neither the 



emitting nor the absorbing gas is expected to share the motion of either star; efforts of this type 
have contradicted each other; and now our FOS4 results create even greater doubts. The observed 
velocities are more likely to represent samples of wind regions near the moving secondary star, 
but in this case one needs to model the complicated varying three-dimensional spatial volumes, 
local absorption may vary, and the FOS4 puzzle still applies. Thus, until more is known about the 
projected velocities, orbit parameters must be based on other considerations. 

The semimajor axis a is non-controversial, 16-19 AU for a total mass in the likely range 
130-220 Mq. Rough limits on the orbital eccentricity e are usually based on the duration of a 
spectroscopic event. Most likely the main parts of an event occur at orbital longitudes within ±90° 
of periastron, i.e., while the two stars are separated less than twice the periastron distance. Since 
the time required to move from -90° to +90° is Aii 80 « (1200 d) x (1 -e) 3 / 2 , eccentricities between 
0.80 and 0.93 seem reasonably consistent with Figure HI Most likely 0.84 < e < 0.91, Aiigo ~ 30 
to 80 days, and the periastron separation is between 1.5 and 2.8 AU. If the first and second A4687 
episodes in 2009 occurred before and after periastron, the most likely time of periastron was close 
tot»0. 

The orbit orientation is more difficult to gu ess. Very probably th e orbit plane is close to the 
Homunculus midplane, with inclination i ~ 45° ([Davidson et al.ll200ll ); any substantially different 
inclination seems unlikely in view of tidal friction near periastron, and would be outside the scope 
of our discussion here. But what orientation does the orbit have within that plane? Most authors 
quote co, the longitude of periastron as defined in textbooks. For rj Car, co = 270° would indicate 
that the secondary star is on the far side of the primary at pe riastron, while co = 1 8 0° represents 



an orb it wh ose major ax i s is p erpendicular to our line of sight. Okazaki et al. (j2008l ). Parkin et al 
(|2009l ) , and iGroh et al.1 (120101 ) all favored co w 240° based on X-rays and other data, but these 



are not strictly independent e stimates, s ince t hey shared a number of insecure assumptions. One 
symptom of uncertainty is that llshibashil (|200ll ) found i ~ 200° by similar reasoning. (The difference 
between 200° and 240° is substantial, because the latter value places the secondary star almost 
behind the primary at periastron, while 200° gives a more "sideways" view.) Ishibashi's analysis 
was far simple r, but t he other, more elaborate analyses required more assumptions. Moreover, 
Kashi Sokerl (|2009bl ) used the X-rays to find an orientation much different from those cited 
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above. Many neglected effects can have major consequences on the wind simulations. For example, 
variable inhomogeneities commonly exist in stellar winds, may be especially crucial for r/ Car's 
X-ray luminosity, and are extremely difficult to model. Therefore u is in fact quite uncertain. 



Another datum may be relevant. iMehner et al.1 (|2010al ) found conspicuous brightness maxima 



in the high-excitation [Ne III] and [Fe III] emission along our line of sight to the star, around 
phase 0.95 or t ~ —110 d in 2003 - i.e., considerably before the event had begun. What is special 
about that part of the orbit? Considering the nature of those emission lines (see Mehner et al.), 
it probably signalled the time when our line of sight passes through the rarified secondary wind - 
which corresponds to the time when the secondary star passed through the projection of our line 
of sight onto the orbit plane. If u = 240° and e < 0.90, this should have occurred at least 160 days 
before periastron. A few calculations show two alternatives for t ~ —110 d: uj ~ 240° and e ~ 0.94, 
which seems excessively eccentric, or uj in the range 220°-232° if 0.85 < e < 0.90. At any rate the 
orientation of the orbit is not well established^ 



Did an eclipse play any role in the spectroscopic event? For the most popular orbit 
parameters (see above), the secondary star should have passed behind the far side of the primary 
wind in the weeks following periastron. The intervening wind might conceivably have been opaque 
in 2003 but not in 2009. For a credible decrease in density between those times, however, the wind 
should not have been very opaque in 2003: the eclipse ingress and egress should not have been 
abrupt. See remarks in lMartin et al,l (|2006al ). 
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Appendix A. Concerning "Phase" and the 5.5-year Period 

We know from experience that varying definitions of "phase" in r] Car's 5.5-year cycle have 
caused confusion. There have been at least three difficulties, some of them matters of principle: 

1. Most authors specified their zero points to to coincide with various observed phenomena, e.g. 
the disappearance of some spectral feature. This policy gives an impression that the definition 
is based on known physical effects, correlated with theory. In fact, the spectroscopic events are 
too poorly understood to identify any observable quantity with a basic physical or geometric 
parameter in the binary system. 



9 Incidentally, the high-excitation precursor peak should recur in April-May 2014, and, if our interpretation is 
valid, it should be broader and less conspicuous than in 2003. Reason: The decreased primary wind should have 
broadened the opening angle of the wind-wind "shock cone." 
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2. This style of reckoning gives a false impression of high precision and future repeatability. In 
fact, successive events are known to differ from each other. Given the continuing changes in 
t] Car (see refs. in Section 1), any observable quantity may shift its phase in the cycle. Terms 
such as "ephemeris" can be very misleading in circumstances like these. 

3. The adopted periods and zero points have varied from paper to paper! We have identified at 
least seven different phase definitions of this type in the literature, with periods spanning a 
range of 10 days and zero points ranging across 20 days. Authors have repeatedly changed 
their definitions. 



Fortunately one simple definition has remained constant, and is greatly preferable in terms 
of logic and procedure. For the rj Car HST Treasury Program Archive ( |http://etacar.umn.edu/[ ), 
"phase" was intentionally defined in a calendar-based rather than a phenomenon-based way. This 
choice was made specifically to avoid the pitfalls listed above. Its period and zero point are 2023.0 
days and MJD 50814.0 = J1998.000 exactly. Phase 2.000 occurred at MJD 50860.0 = J2009.077. 
The integer quantities help to minimize calculative errors, and, more important, they discourage 
any impression that t = represents some critical time with respect to physics. (Periastron is 
thought to occur fairly close to t = 0, but this categorically plays no role in the definition.) 

This system has been in use since 2003 without alteration; its ado pted period of 2023.0 



days continues to agree wit h the best measurements within less than la (jDamineli et al.ll2008bl ; 



Fernandez-Lajus et al.ll2010l ): and it is used in the Treasury Program archive, which is the largest 
easily-accessible source of data on r] Car (reduced data of all HST STIS, VLT UVES, and Gemini 
GMOS observations are available). Hence there is no reason to substitute any later, arbitrary sys- 
tem. In order to minimize present and future confusion, and for the other reasons noted above, it 
is the obvious standard. 

Timing measurements are valuable because they may indicate changes between successive 
events. One particular operational procedure achieves the best precision and reproducibility, as 
follows. Various observable quantitities - e.g., photometry at most wavelengths - attain brief 
maxima Qmax at a particular stage in a spectroscopic event, and then briskly fall to their local 
minima Qmin (or vice- versa). The times of maximum and minimum are imprecise, but for some 
observables the time of midpoint, when Q is the average of maximum and minimum, is more precise 
than anything else in the spectroscopic event! This is true mainly because the descent is very rapid 
at that time. If Q(t) is fairly smooth and if enough data points are available, then the most robust 
measurement protocol is as follows. 



1. Estimate the value of Q mSiX but ignore its time t(Q m ax) which may be ill-defined. A local 
quadratic fit may be appropriate, but fortunately there is no need for very high precision in 
Qmax (see below). 

2. Do the same for Q m i n , and calculate the midpoint Q m = (Qmin + Qmax)/2. 
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3. Then estimate the midpoint time t m = t(Q m ), based only on the data points that are nearest 
to it. If many data points are available, for instance, one might use only those in the range 
Q m ± AQ where AQ « (Q 

max — Qmin) / 4. In favorable circumstances either a linear or a cubic 
fit gives practically the same result. 



In effect this is a specialized form of interpolation. Each step employs only a local subset of the 
data, so local irregularities in behavior around the maximum and minimum do not affect t m much. 
In the case of rj Car's events, t m is only weakly sensitive to the estimated values Q max and Q mm . 
This is true because, for most observables, the curve Q(t) is approximately antisymmetric near the 
midpoint and the rate of descent is quite rapid there. Note that theoretical significance plays no 
role in the e mpirical reasoning. For an especially successful example, r] Car's 'J' magnitudes in 2003 
reported by lWhitelock et al.1 (120041 ). the rms formal error in t m is only ±0.4 day. We employed the 
same method for some details in Section 14.11 above. 



Fernandez-Lajus et al.l (|2010l ) used a method that is more general but is not as well suited to rj 
Car in particular. In effect, their procedure mixes data in a broader time interval, thereby allowing 



a stronger dependence on behavior details around Qmax 

and Qmin. iDamineli et al.1 l2008bl defined 
a phase zero point based on an extrapolation, which of course is inherently far less robust than 
interpolation. 



Fernandez-Lajus et al.l (|2009l . |2010| ) report high-quality ground-based pho tometry sinc e 2003 , 
but it is extremely unfort unate that the earlier JHK program described by iFeast et al.l (|200ll ) 
and IWhitelock et al.l (J2004J) was forced to close before the 2009.0 event. Their results give precise 
estimates of a time related to the 2003.5 event, and their earlier data give fairly good results for the 
1998.0 event (see below). Here we estimate the 2003-2009 inter-event time interval from our HST 
UV photometry (Section [3]), because these data represent the stellar wind with little contamination 
by ejecta at r > 0.15". Unfortunately the temporal sampling is sparse compared to ground-based 
photometry. 

Table H] lists the results for t m in three events and five wavebands. Here we use JHK photometry 
by Feast et al. and Whitelock et al. cited above, and our UV photometry. JHK wavelengths 
are pertinent because light reflected by dust in the ejecta is less than for UBVRI. The quoted 
uncertainties are semi-formal rms standard errors based on observed scatter or deviation of the 
data points and other details; these appear to be realistic, judging from the scatter in the three 
periods based independently on J, H, and K. Note that t m for a particular event may depend on 
wavelength; it is the difference between events that we are concerned with here. The JHK data 
indicate an interval 202 1.5 ± 0.9 d between t he 1998.0 and 2003.5 events; this compares with 
2022.7 ± 1.3 d found by IDamineli et al.l (2008b) with less suitable methods and a much larger set 
of data, and 2022.8 ± 0.5 d found bv lFernandez-Laius et al.l (|2010l ) from ground-based photometry 
in 2003-2009. The HST 330 nm value is obviously consistent with these, but the 250 nm result 
appears discordant. This may be a real effect; the 250 nm wavelength region physically differs from 
the others because it represents the Fe II forest (Section [3]). Unfortunately there is no assurance 
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that UV photometry will be possible during the next, 2014.6 event. 
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Fig. 1. — HST UV photometry. The upper panel shows the light curve from 1998 to 2011 (boxes: 
ACS/HRC, stars: WFPC2, crosses: STIS/CCD, values for F330W are shifted by -0.4 mag). The 
two lower panels compare F250W and F330W brightness variations during the 2009.0 spectroscopic 
event (filled circles and dotted line) with the 2003.5 spectroscopic event (open circles). Day "0" 
corresponds to MJD 52837 for the 2003.5 event and MJD 54860 for the 2009.0 event. The black mark 



i n the lower left panel indicates the date of the V band minimum observed by lFernandez-Lajus et al 

(boid ). 
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Fig. 2. — Time sequence of the He II A4687 emission during the 2009 event in Gemini GMOS 
observations. Continuum was normalized to unity at A4740 and is indicated with horizontal dotted 
lines on the right side. Offset between tracings is 0.1. Number of days before (— ) and after (+) 
MJD 54860 are indicated next to each spectrum. The dotted vertical line indicates the position of 
He II A4687 at zero radial velocity. 
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Fig. 3. — Equivalent width and radial velocity measurements of the He II A4687 emission on the 
star during the 2009 event in Gemini GMOS data. The dotted vertical line indicates the time of 
maximum negat ive radial velocity which o ccurs at the flux-decline-midpoint. This is also true for 
previous events (jSteiner &: Daminelil 12004). 
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Fig. 4.— TOP: X-ray light curve during the 1998, 2003.5, and 2009 events (from 
|http: / / asd.gsfc.nasa.gov/Michael.Corcoran/eta_car/etacar_rxteJ ightcur ve/index.htmll ) . BOT- 
TOM: Equivalent width of the He II A4687 emission during the 2009 even t with Gemini GMOS 
(black squares) and selected data points (black triangles) from lTeodoro et al.l (120111 ) which constrain 
the timing of the second episode during the 2009 event further. The red squares are measurements 
of HST STIS data during the 2003.5 event. No significant He II emi ssion was found in 2003 at 
t ~ +20 d in STIS data and measurements bv lSteiner Daminelill2004j during the 1992.5 and 1998 
events (green triangles; these values are scaled by a factor of 3, see Footnote [5]) ind icate a delayed 



second episode. This casts doubt on the timing of the second episode found by iTeodoro et al 
(|201ll ). Values for older events are shifted by multiples of 2023 days. 
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Fig. 5. — Equivalent width and radial velocity of He II A4687 at F0S4 during the 2009 event 
observed with Gemini GMOS (red squares). Dashed curves are measurements on the star in direct 
view. The time-delay between spectra of the star in direct view and spectra at FOS4 is about 
18 days. Re-measurements of VLT UVES data during the 2003.5 event are consistent with this 
finding (blue triangles, values are shifted by 2023 days). 
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Fig. 6. — Equivalent width and radial velocity of the He I A4714 emission (filled squares) on the 
star during the 2009 event observed with Gemini GMOS. The radial velocity of the absorption line 
is also shown (filled triangles). The crosses show the values for He II A4687 emission; the correlation 
is obvious. 



-40- 



Wavelength 

4090 4100 4110 4090 4100 4110 4090 4100 4110 4090 4100 4110 

4.5 n 1 n pt 1 ■— i 1 rn 1 rn 




-1000 -500 -1000 -500 -1000 -500 -1000 -500 

Radial Velocity (km s" 1 ) 



Fig. 7. — H(5 in tracings across the SE lobe in observations at t = —353 d (2008-02-11, phase 
1.83), t = -82 d (2008-11-08, phase 1.96), t = +21 d (2009-02-19, phase 2.01), and t = +344 d 
(2010-01-08, phase 2.17). The key indicates the positional offset along the spectrograph slit, SE 
from the central source. Spectra are shifted to compensate for their A^ "moving mirror" redshifts. 
The red curve (4'.'5) corresponds to a spectrum which is close to position FOS4. 
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Fig. 8. — Tracings across the SE lobe showing He I A4714 emission and absorption at the same 
times as Figure Velocities are corrected for the "moving mirror" redshift AV. The red curve 
(4'.'5) refers to position FOS4. Note that the large blueshift of the He I emission and absorption 
lines during the events, discussed in Section 5, cannot be observed in this Figure since it occurs 
between t « — 70 d and t ~ d not sampled here . The feature at the left edge is [Fe III] A4703, 
which disappears during each spectroscopic event (IMehner et alJl2010al ). 
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Fig. 9. — Absorption components of Fe II and similar ions at FOS4 deepened for a few months 
during the 2009 event (solid red curve). They can be identified with absorption components at — 
400-450 km s _1 of Fe II lines (solid black bars) and Cr II, Mg I, and Ti II lines (dotted bars). 
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Table 1. HST Photometry Results' 









STM Afi b 


MJD 


Year 


/in — 11 pro . rrn — 2 o — 1 A— 1^ 


(mag) 






WFPC2 PC/F255W Filter 




54717.1 


2008.69 


1.074 


6.323 


54787.8 


2008.88 


1.126±0.013 


6.271±0.013 


54807.4 


2008.93 


1.099±0.035 


6.299±0.035 


54834.0 


2009.01 


1.189±0.004 


6.212±0.004 


54838.4 


2009.02 


0.967±0.014 


6.437±0.016 


54841.5 


2009.03 


0.806±0.006 


6.635±0.008 


54845.5 


2009.04 


0.598±0.017 


6.959±0.031 


54854.5 


2009.06 


0.451±0.003 


7.265±0.008 


54872.4 


2009.11 


0.437±0.019 


7.298±0.048 


54949.2 


2009.32 


0.937±0.011 


6.471±0.013 






WFPC2 PC/F336W Filter 




54717.1 


2008.69 


1.108±0.035 


6.289±0.034 


54787.8 


2008.88 


1.114±0.041 


6.283±0.040 


54807.4 


2008.93 


1.142±0.031 


6.256±0.030 


54834.0 


2009.01 


1.351±0.005 


6.073±0.004 


54838.4 


2009.02 


1.399±0.037 


6.036±0.029 


54841.5 


2009.03 


1.388±0.023 


6.045±0.018 


54845.5 


2009.04 


1.300±0.007 


6.116±0.006 


54854.5 


2009.06 


1.004±0.020 


6.396±0.022 


54872.4 


2009.11 


1.008±0.027 


6.392±0.029 


54949.2 


2009.32 


1.599±0.009 


5.890±0.006 




STIS Synthctic/F250W Filter 




55062.0 


2009.63 




6.12 


55258.7 


2010.17 




5.89 


55428.3 


2010.63 




5.65 




STIS Synthetic/F330W Filter 




55062.0 


2009.63 




5.63 


55258.7 


2010.17 




5.55 


55428.3 


2010.63 




5.35 



"Earlier results c an b e found in llMartin fc Koppelmanl 120041; 
iMartin et aTll2006bl. bold) . 

b The STM AG photometric system is calibrated for direct com- 
parison of fluxes with similar filters in different instruments. 



Table 2. Gemini GMOS Journal* 











Star c 


FOS4 d 


Name b 


Date 


MJD 


Phase Ccnwavc 


Exp. timc e 


Exp. timc e 








(A) 


(s) 


(s) 



gH45 


Jun 16, 2007 


54267.1 


1.707 


4300 


10 


- 


gH45 


Jun 17, 2007 


54269.0 


1.708 


4300 


10 


- 


gH49 


Jun 30, 2007 


54281.0 


1.714 


5600 


40 


220 


gill 


Feb 11, 2008 


54507.4 


1.826 


4300 


77 


377 


gill 


Feb 13, 2008 


54509.2 


1.827 


5600 


4 


260 


gI50 


Jul 5, 2008 


54652.0 


1.897 


5600 


40 


260 


gI54 


Jul 17, 2008 


54665.0 


1.904 


4300 


10 


300 


gI85 


Nov 8, 2008 


54778.3 


1.960 


5200 


40 


453 


gI85 


Nov 8, 2008 


54778.3 


1.960 


4300 


77 


377 


gI90 


Nov 27, 2008 


54797.3 


1.969 


5200 


40 


450 


gI90 


Nov 27, 2008 


54797.3 


1.969 


4300 


77 


377 


gI96 


Dec 18, 2008 


54818.3 


1.979 


5200 


40 


450 


gI96 


Dec 18, 2008 


54818.4 


1.979 


4300 


7 


377 


gI98 


Dec 25, 2008 


54825.3 


1.983 


5200 


40 


450 


gI98 


Dec 25, 2008 


54825.4 


1.983 


4300 


77 


377 


gI99 


Dec 31, 2008 


54831.3 


1.986 


5200 


40 


450 


gI99 


Dec 31, 2008 


54831.4 


1.986 


4300 


10 


70 


gJOl 


Jan 4, 2009 


54835.3 


1.988 


5200 


45 


450 


gJOl 


Jan 4, 2009 


54835.3 


1.988 


4300 


55 


377 


gJ02 


Jan 9, 2009 


54840.2 


1.990 


5200 


45 


450 


gJ02 


Jan 9, 2009 


54840.2 


1.990 


4300 


33 


377 


gJ03 


Jan 12, 2009 


54843.3 


1.992 


5200 


45 


450 


gJ03 


Jan 12, 2009 


54843.3 


1.992 


4300 


6 


377 


gJ04 


Jan 15, 2009 


54846.2 


1.993 


5200 


15 


450 


gJ04 


Jan 15, 2009 


54846.2 


1.993 


4300 


33 


377 


gJ05 


Jan 21, 2009 


54852.3 


1.996 


5200 


4 


450 


gJ05 


Jan 21, 2009 


54852.3 


1.996 


4300 


6 


300 


gJ06 


Jan 24, 2009 


54855.3 


1.998 


5200 


4 


260 


gJ06 


Jan 24, 2009 


54855.4 


1.998 


4300 


6 


300 


gJ07 


Jan 29, 2009 


54860.4 


2.000 


5200 


9 


450 


gJ07 


Jan 29, 2009 


54860.4 


2.000 


4300 


33 


377 


gJ09 


Feb 5, 2009 


54867.2 


2.004 


5200 


15 


450 


gJ09 


Feb 5, 2009 


54867.3 


2.004 


4300 


6 


300 


gJ13 


Feb 19, 2009 


54881.2 


2.010 


5200 


4 


450 


gJ13 


Feb 19, 2009 


54881.3 


2.011 


4300 


6 


300 


gJ20 


Mar 17, 2009 


54907.3 


2.023 


5200 


4 




gJ20 


Mar 17, 2009 


54907.3 


2.023 


4300 


6 




gJ32 


Apr 28, 2009 


54949.1 


2.044 


5200 


4 


260 


gJ32 


Apr 28, 2009 


54949.1 


2.044 


4300 


6 


300 


gJ56 


Jul 23, 2009 


55036.0 


2.087 


5200 


4 


260 


gJ56 


Jul 24, 2009 


55036.0 


2.087 


4300 


77 


300 


gK02 


Jan 8, 2010 


55204.32 


2.170 


5200 


4 




gK02 


Jan 8, 2010 


55204.34 


2.170 


4300 


4 




gK05 


Jan 20, 2010 


55216.29 


2.176 


5200 




260 
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Table 2 — Continued 



Name b Date MJD Phase Ccnwave 

(A) 


Star c 
Exp. timc c 
(?) 


FOS4 d 
Exp. timc c 





a PA = 160°, Grating: B1200.G5321. 
b As listed in http://etacar.umn.edu/ 

c Include slit positions less than ±0"375 from the star at A4687. 

d Include slit positions -1 to -3" from the star at A4687. 

c Several exposures were taken on each data, combined exposure times are 
listed. 



Table 3. Equivalent width and radial velocity of He II A4687 in Gemini GMOS data on the star 



Name* 


Date 


MJD 


Phase 


EW 


rad 




CUT) 






(A) 


(km s" 1 ) 


gH45 


2007 Jun 16 


54268.0 


1.707 


-0.25 ± 0.02 


15.09 ± 6.11 


gH49 


2007 Jun 30 


54281.0 


1.714 


-0.19 ± 0.01 


-39.60 ± 46.88 


gill 


2008 Feb 11 


54507.4 


1.826 


-0.34 ± 0.03 


-35.39 ± 54.12 


gI50 


2008 Jul 05 


54652.0 


1.897 


0.01 ± 0.06 


-86.80 ± 10.60 


gI54 


2008 Jul 17 


54665.0 


1.904 


-0.25 ± 0.01 


-50.83 ± 40.44 


gI85 


2008 Nov 08 


54778.3 


1.960 


0.82 ± 0.30 


-50.46 ± 12.29 


gI90 


2008 Nov 27 


54797.3 


1.969 


0.60 ± 0.23 


-151.68 ± 2.73 


gI96 


2008 Dec 18 


54818.3 


1.979 


2.16 ± 0.23 


-209.59 ± 21.47 


gI98 


2008 Dec 25 


54825.3 


1.983 


2.09 ± 0.15 


-170.93 ± 21.29 


gI99 


2008 Dec 31 


54831.3 


1.986 


2.16 ± 0.17 


-279.30 ± 7.06 


gJOl 


2009 Jan 04 


54835.3 


1.988 


2.76 ± 0.19 


-311.94 ± 13.91 


gJ02 


2009 Jan 09 


54840.2 


1.990 


2.40 ± 0.21 


-328.94 ± 4.58 


gJ03 


2009 Jan 12 


54843.3 


1.992 


1.21 ± 0.23 


-418.11 ± 12.49 


gJ04 


2009 Jan 15 


54846.2 


1.993 


0.24 ± 0.22 




gJ05 


2009 Jan 21 


54852.3 


1.996 


0.00 ± 0.18 


-118.72 ± 30.72 


gJ06 


2009 Jan 24 


54855.3 


1.998 


0.20 ± 0.17 


-114.85 ± 49.71 


gJ07 


2009 Jan 29 


54860.4 


2.000 


0.34 ± 0.23 


-99.95 ± 59.14 


gJ09 


2009 Feb 05 


54867.2 


2.004 


0.61 ± 0.18 


-95.55 ± 17.06 


gJ13 


2009 Feb 19 


54881.2 


2.011 


1.10 ± 0.19 


3.08 ± 3.97 


gJ20 


2009 Mar 17 


54907.3 


2.023 


0.22 ± 0.25 


64.57 ± 59.78 


gJ32 


2009 Apr 28 


54949.1 


2.044 


0.14 ± 0.26 


9.28 ± 33.48 


gJ56 


2009 Jul 24 


55036.0 


2.087 


0.06 ± 0.30 


27.39 ± 19.36 


gK02 


2010 Jan 08 


55204.3 


2.170 


0.08 ± 0.19 


-56.88 ± 7.80 


a As listed on the Eta Carinae Treasury Project site at|http:/ /etacar.umn.edu/| 
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Table 4. Photometric t p f lot MJD's for rj Car's 1998.0, 2003.5, and 2009.0 spectroscopic events 





1998.0 


2003.5 


2009.0 


Differenz 


Waveband 


tphot 


tphot 


tphot 




K a 


50802.6 ± 1.8 


52824.7 ± 0.5 




2022.1 ± 1.9 d 


H a 


50803.3 ± 1.3 


52825.9 ± 0.4 




2022.6 ± 1.4 d 


J a 


50807.0 ± 1.5 


52826.8 ± 0.4 




2019.8 ± 1.6 d 


330 nm 




52827.8 ± 1.7 


54850.0 ± 1.7 


2022.2 ± 2.4 d 


250 nm 




52826.1 ±1.1 


54842.7 ± 0.6 


2016.6 ± 1.3 d 



HFeast et all l|200ll ); IWhitelock et all l|2004h . 



